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Abstract 

Within the framework of variational modelling we derive a two-phase moving boundary 
problem that describes the motion of a semipermeable membrane separating two viscous 
liquids in a fixed container. The model includes the effects of osmotic pressure and surface 
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1 Introduction 

This paper is devoted to a two-phase moving boundary problem describing osmotic swelling of 
a closed membrane in a viscous liquid. 

Let C C M. N (N > 2) be a bounded connected open set with smooth boundary representing a 
fixed region filled with an incompressible viscous liquid that moves according to the velocity field 
u = u(t,x). Inside the liquid there is a closed connected semipermeable membrane T(f) C C 
enclosing an open set 17+(t) and separating it from the outer phase 12_(i) := C\Ct + (t). In both 
phases a certain amount of a solute is dissolved. Its scalar concentration c = c(t , x ) evolves 
by convection along u and diffusion through the liquid. It may be discontinuous across T(t). 
Both the diffusivities and the viscosities are assumed to be constant and positive but possibly 
different in 17+ and 17_, respectively. 
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The membrane is permeable for the liquid but impermeable for the solute. Its deformation 
and movement are governed by surface tension forces, osmotic pressure, and the fluid motion. 
Based on these assumptions the following moving boundary problem can be derived using the 
approach of variational modelling, see Section 2: 


-i^±Am ± + V(g± + c±) 

= 0 

in n±(t), 

t > 0, 

divM-t 

= 0 

in Q±(t), 

t > 0, 

{T(u,q + c)]n 

= Hn 

on T(t), 

t > 0, 

N 

= 0 

on r(t), 

t > 0, 

U- 

= 0 

on dC, 

t > 0, 

d t c± — k±Ac± + Vc± • u± 

= 0 

in 

t > 0, 

n±d n c± + c±([c] + H) 

= 0 

on r(t), 

t > 0, 

d n c~ 

= 0 

on dC, 

t > 0, 

V n 

= H + [c| + u ■ n 

on T(t), 

t > 0, 

F(0) 

= r°, 



c±(0) 


in n±(0), 

> 


where we used the notation u± := u |n ± , c± := c|n ± - The brackets [•] indicate the jump of a 
quantity across T(i), i.e. 


MV)](aO := I™ w(t,y) - lim w{t,y) 

y(z£l+(t),y—tx y£il—\t),y—>x 


for x £ T(f) and w : Uf2_(f) — > K. Further, H = H(t,x ) denotes the (N — 1) - fold mean 

curvature of the closed compact hypersurface T(t) = df,l+(t) at the point x £ T(f), oriented 
in the way that spheres have negative curvature, while V n is the normal velocity of the family 
(r(t)} w.r.t the unit normal field n = n(t) of T(t) pointing outward The operator d n 

takes the directional derivative of a sufficiently regular function w.r.t the normal field n(t). If 
no confusion seems likely, we use the same symbol d n to denote the derivative in the direction 
normal to dC and exterior to C as well. The symbol q = q(t,x) stands for the hydrodynamic 
pressure and 

r±{u±,q ) := ^ ± e(M±) - q± Id := ^ ± (Vu± + (Vu ± ) T ) - q± Id, q± := g|n ± , 


is the hydrodynamic stress tensor. Note that the initial velocity u( 0) is uniquely determined 
by c(0) and T(0) as we shall discuss later in some detail, cf. Section [7T1 

System (ED is written in dimensionless form. The given positive constants k± and v± carry 
information about physical parameters such as diffusivity of the solute, viscosity of the liquid 
in both phases and permeability of the membrane to solvent. For later use we introduce the 
piecewise constant functions 


v(x, t) 


v+ x £ £l+(t), 
v- x £ fl-(f), 


k(x, t) 


K + X £ 
K- X £ 


In the corresponding one-phase situation, a detailed derivation of the model within the frame¬ 
work of variational modelling has been given in [10] - The two-phase problem is obtained in 
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a parallel fashion. Therefore we restrict ourselves here to a brief recapitulation of the chosen 
setup which will be given in Section [2] 

The paper m also contains a short-time existence result for classical solutions for the 
one-phase problem. For the simpler limit problem in which the membrane moves through an 
immobile liquid the existence of classical solutions for a short time has been established in 

and the paper [5] deals with a stability analysis of its equilibria. For different modelling 
approaches (excluding fluid motion) as well as analytic results in even more special situations 
such as radial symmetry we refer to nmnniip and the references given in uni. 

In this paper we focus on the equilibria of system m and their stability properties. These 
equilibria form a finite dimensional submanifold of the phase space. Our main result states that 
this manifold is locally exponentially attractive, i.e the system is normally stable in the sense 
of [14] . 

While the main line of the proof is parallel to the one in !9], we have to deal with the 
additional difficulty of handling the nonlocal solution operator of the two-phase Stokes system 
(flll i - (0)5. In particular, the results of [!] that are a crucial ingredient of the stability 
analysis in BUM are no longer directly applicable. Additionally, one has to discuss the full 
two-phase Stokes system with respect to well-posedness and regularity. As these results do not 
seem to be readily and explicitly available in the literature, we include a proof of them in an 
appendix. 

The present paper is organized as follows. In Section [2] we explain briefly how our model 
can be derived within the framework of variational modelling. Section [3] identifies the equilibria 
of system m- In Section 2] we transform the problem to a fixed reference domain, determine 
the linearization of the transformed problem around an arbitrary fixed equilibrium and analyse 
some spectral properties of the corresponding linear operator. In this section we also give 
a precise formulation of our main result ("Theorem 14.21) . which is proved in Section [5] The 
appendix (Section 0 contains a detailed discussion of the full two-phase Stokes system (17.11) 
and some abstract facts that are helpful for the spectral analysis (Lemma 15.11 Corollary 15.21 
Lemma fOll . 

2 Modelling 

We use the same modelling approach as in [TO] and derive our model from the following building 
blocks: 

i) We consider paths in a state manifold Z consisting of pairs (fl + , c) of a simply connected 
domain 0_|_ satisfying C C and a nonnegative solute concentration c : C —> R that may 
be discontinuous across <9fl + . The domain and the container C uniquely determine 
0 _. 

ii) On Z we define the energy functional 

£(fl + ,c ):=7 f clwcdx + a\d£l + \ ( 2 - 1 ) 

Jc 

with positive constants a and 7 , cf. . 10 ) for a more detailed discussion of their physical 
meaning. This choice includes diffusion of the solute and surface tension of the membrane 
as driving mechanisms of the evolution. 
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iii) The processes that dissipate energy are solvent motion, solute flux, and passage of solvent 
through the membrane. Taking into account incompressibility of the solvent and mass 
conservation of the solute, these processes can be represented by triples 

{(it, /, V n ) | divu = 0 in C, [wfl = 0 on <90+, f±-n = c±V n on <90+} 
which we collect in the process space Vm,, <■)■ 

iv) The dissipation functional is defined on V(q, + )C ) and given by 

+ ,c)(u,f,V„):=l- f —-—— dx + \- f u 2 \e(u)\ 2 dx + ^ [ (u ■ n - V n ) 2 da, 

1 Jc c z JG 1 J8Q.+ 

where v.j = in 0 + (j = 1 , 2 ) and 1/3 are positive constants related to mobility of the 
solute, viscosities of the solvent in both phases and to the membrane’s permeability to 
the solvent, cf. again m- 

v) Observe that the elements of the tangent spaces C \Z of the state manifold Z can 
be represented by pairs (V n ,c), where V n : 90 + —> K is a normal velocity and c is a 
concentration change. Since mass conservation of the solute is expressed by the relation 
c + div/ = 0 , it seems natural to define the process map II(q + jC ) : 7 :> (n + , c ) —t T(n + ,c)Z 
defined by 

n (Q +! c)(u,f,V n ) = 04, -div/). 

The model m is now determined by the dynamical system on Z 

z = Il z w*, z=(Q + ,c), (2.2) 

where w* is the solution to the minimization problem 

z {w) + £'(z)\R z w] —*• min, ui€V z , (2.3) 

cf. [ID], and by an appropriate scaling. 


3 Equilibria 

Observe that by construction the system (EJ is a gradient flow w.r.t the functional 

£(fi+,c)= / clncda: + |r| 

Jc 

(cf. pH] Section 2 for a more detailed discussion of this fact). Hence, the functional £ is a 
Ljapunov function for the system ED. Indeed, assuming smoothness and strict positivity of 
concentrations, integration by parts yields 

| £(n+(i), c(i)) = - - J m) 1 + m? -\f c K«(i))| 2 - (3.1) 
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Let (u, q, c, I2 + ) be an equilibrium solution to m (i .e. (u, q, c,fl+) is constant in time, T is a 
closed connected hypersurface, and u,q,c are continuously differentiable away from T). Since 


m vanishes at equilibria, Korn’s inequality implies that u = 0. Moreover, c must be constant 


in both phases and [c] = —H. Thus, also H is constant, so that T is a sphere. The first 
equation in O) implies then that q is constant in both phases, and from the third equation 
one concludes that [g|n = 0 on T. Summarizing: 

Lemma 3.1 A tupel (zt, g, c, f2+) is an equilibrium solution to 0) iff 12-|- is a ball of some 
radius R, f2_|_ C C, [c] = ( N — 1 )/R, u = 0 and q is constant in C. 

4 Linearization at an equilibrium 

We fix now a single equilibrium (0 ,q,c,D + ) and assume w.l.o.g. that D + = B(0,1) and [c] = 
N — 1 =: m. We further define S := c®(0, 1), D_ := C \ D + and keep these notations fixed 
hereafter. 

In order to solve system (HD we are going to consider a set of transformed equations given 
over D± as fixed reference domains. The unknown family of surfaces |F(t)} will be described 
by a signed distance function with respect to the unit sphere. The ansatz is standard and has 
already been used in [9] in an identical way. 

The mapping 


X : S x (—1,1) —> M . n , (x, s) (1 + s)x, 


is a smooth diffeomorphism onto its range. Fix 0 < a < 1 small enough that D C C, where 
D := range(X|s X ( a a )). As it is convenient, we decompose the inverse of X := A'|g X (_ aja ) into 
X~ x = (P, A) : D —> S x (—a, a), where P is the metric projection onto S, and A is the signed 
distance function with respect to S, i.e. P(x) = x/|x|, A(x) = |x| — 1. Let a £ (0, a/4) and 


Ad := {cr e C 1 (S); ||cr|| C (S) < a}- 


It is well-known that, given a £ Ad, the mapping 9 C7 {x) \= (1 + a{x))x is a diffeomorphism 
mapping S onto S a := 0 a [S] ■ We extend this diffeomorphism to the whole of R N : Let \ G 
C°°(]R, [0,1]) satisfy x\[-d,d] = 1 , xl(-oo,-3o] = xl[ 35 ,oo) = 0, Hx'lloo < 1/a- Then the mapping 



X{P(y),A{y) +x(A(y)) -a(P(y)), if y £ D 


(4.1) 


if y (£ D 


y, 


(a £ Ad), again denoted by 9 a , is an appropriate extension, the so-called Hanzawa diffeomor¬ 
phism. We have 9 a £ Diff(R Ar , 'R N ). Moreover, 9 a = id outside D , in particular in a sufficiently 
small open neighborhood of dC. Moreover, denoting by D a , + the domain enclosed by S a and 
letting Da,- := C \ D a+1 we have that 


9a\ D± GDiff(£>±,£> ff ,±), 


cr £ Ad, and dD G ,+ = S a , dD a - = S a U dC. Finally note that the surface S a is the zero level 
set of the function ip a defined by 


tfia{x) = A(x) - cr(P(x)), 
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x £ D, a £ Ad, i.e. 5 CT = </5 cr 1 [{0}]. For later use we set 


L a (x) \= \S7ip a \(9 a (x)). 


It can be shown that L a > 0 on S for all a £ Ad. 

Given a £ Ad, let 9*, 91 denote the pull-back and push-forward operators induced by 9 a , i.e. 
9* f = fo9 a , 91 g = go9~ x . If the functions 6,p are time dependent, i.e. b = b(t,x ), p = p(t,x), 
we define [9* b\(t,x) := [9b(t , -)](a;), analogue for 91. 


Using this notation, for p : J C [0, oo) 
introduce the transformed operators 


n{p) 
H(p ) 
A±(p)w± 
B±(p)w± 
IC±{p)w± 


Ad n C 2 ( S ) and sufficiently smooth w± £ l c± we 


e *P n ls P ]\ 

9 *p H [S P ]5 

e;(A {0Sw ± )y, 
0*p(y(9lw±)\ Sp ) ■ n{p); 
o;(v(9iw ± )). 


Letting p± := c± o 9 p , '■= c±.o o 9 Po , instead of (11.11) . we study the following problem on 

D± as fixed reference domains: 


d t g± - n±A(p)g± +/C±(p)/i± ■ s±(p) + R±(p,p) 
k±B(p)p± + p±(M +H(p)) 
dnP— 

d t p - L(p)[H(p) + [p] + s{p) ■ n] 

P±(0) 


P( 0) = 


0 

0 

0 

0 

P±, o 

Po, 


in D±, 
on dD + , 
on <9C, 
on <9.D + , 
in L>±, 


(4.2) 


where s(p) := 9*u , 


+ Vg± = 0 

div u± = 0 

\r(u,q)\n = Hn 

N = o 

U- = 0 


in D p ±, t > 0, 
in D Pi ±, t > 0, 
on Ap, t > 0, 
on S p , t > 0, 
on <9C, t > 0, 


(4.3) 


and s±(p) := s(p) |z? ± ■ The terms 7?.± arise from the transformation of the time derivative ( p±)t 
and are determined by 


H±(w±,a)(y) = r 0 (L a [H(a) + [w]| + s(p) ■ n(p)],B li (a)w±)(y ), 
where w± £ C 1 (D±), a £ Ad and 


V G D±, 


r 0 (h,k)(y) := 


X(A (y)) ■ h(P(y)) ■ k{y ), if y £ D 
0, if y £C\D, 


(4.4) 


B p (a)w±(y) = 9* S7(91w±)(y) ■ (n s o P)(y), y £ D± 

(ns being the exterior unit normal field of S). The explicit calculation 
of 1Z± is straightforward, cf. again [5]. 
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Linearization of (14.21) around the equilibrium (p,±,p ) = (c±,0) yields the following system 
for the shifted variable /i — c, denoted again by p: 


d t p± - k ± Ap± 

= F±(P±, P ) 

in D±, 

K±d n p± + c±([/z] + (As + m)p) 

= G±(p±,p) 

on 5, 

dnP— 

= 0 

on 9G, 

dtp - [(As + m)p + [/z] + s'(0)p ■ n] 

= H(p ±1 p) 

on 5, 

p±( 0 ) 

= P±, 0 

in D±, 

P(0) 

= Po, 



with suitable nonlinear remainders F, G, H that act smoothly between the function spaces we 
are going to use, cf. Lemma 4.2 in [8] and Corollarv l7.2l in the present paper. By construction, 
they satisfy 

F±( 0) = G±(0) = H( 0) = 0, i^(0) = G' ± (0) = H'( 0) = 0. 

By A s we denote the Laplace-Beltrami operator of the unit sphere. After some algebra, letting 
a± := k±/c±, A := As + m and 

Gi(p±,p) = G + (p±,p)/c + ; 

G 2 {p±? p) G + (p±, p)/c + G_ {p ±, p)/c —, 

G 3 (p±,p) = H(p±,p) +Gi(p±,p), 


we get 


d t p± - n±Ap± 
a+d n p+ + \pj + A p 
[ctdnpj 
dnp— 

d t p + a+d n p+ - s'(0)p ■ n 
P±( 0) 
P( 0) 


= F±(P±,P) 

in D±, 

= Gi(p±,p) 

on 5, 

= G 2 (p±,p ) 

on 5, 

= 0 

on dC , 

= G 3 (p±,p) 

on S , 

= P±, 0 

in D±, 

= Po- 



(4.6) 


We close this section by defining the abstract setting for our analysis and giving a precise 
statement of our main result. Let p > N + 2. For s > 0 and a Banach space Y, M £ 
{D±, S, [0, T], [0, 00)} (T > 0) we denote by W£(M,Y) the L p -based Sobolev space of order 
s. In particular, if s £ N, this fractional-order Sobolev space coincides with the Besov space 
Bp p (M,Y) (cf. [T7j). For the sake of brevity we write W p (M) := W p (M,R) and introduce the 
notations W p (D±) := W p (D + ) x Wg(D-), {p±,p) := (p + ,p-,p). Let 


El := {(p±,p) £ Wp(D±) x W*- 1/P (S); d n p _ = 0 on dG}; 

E := {{p±,p) £ Wp~ 2/p (D±) x W*~ 3/P (S)- d n p _ = 0 on c>G}; 

E 0 := LP(D±) x W 3 ~ 1/p {S) 


and for an interval J C [0, 00 ) 

E(J) := L p (J,Ei) n (Wt(J,IS(D ± )) x W^-^^L^S))). 
We further define spaces of exponentially decaying functions 

E(tf) := {(^±, a) £ E(R+); e 5t (£±, a) £ E(R+)} 


7 






(S > 0 ), equipped with the norm ||(^±, cr) ||e( 5 ) := H el 5 t (£±; ct )IIe(r+)- and recall the standard 
embedding result 

(4.7) 

We formally introduce the operators L, K and B by their action as follows: 

L{p±,p) := (K±Ap±,-a+d n p+), K(p±,p) := ( 0 , 0 , s'(0)p|s • n), 


B(p±,p) ■■= ( a+d n p+ + M + A p, [ ad n pj), 
where u = u(p) = s'(0)p (and a suitable p ) solve 


- Vp± 

= 0 

in D±, 

divu± 

= 0 

in D±, 

[u,p)jn 

= Apn 

on S, 

N 

= 0 

on S, 

U- 

= 0 

on dC. 


(4.8) 


Then, with F := (F±, G 3 ), G := (Gi, G 2 ), (14.61) can be written as an abstract evolution problem 
d t p - {L + K){p) = F(p), Bp = G{p), p(0) = po := (p±,o,Po), p:={p±,p). (4.9) 

Solutions to (14.91) are paths in the manifold 

A4 := {p G E 5 Bp = G(p)}. 

They are supposed to possess the following regularity: 

Definition 4.1 A global strong solution of the evolution problem m is a solution p = 
{p±,p) : [0,oo) —>• E such that 


p\[0,T\ eE([0,T]) V T > 0. 

Observe that the set of equilibria of (14.91) is 

£'■={£& M; —(L + K, B)(e) = (F, G)(e)}, 

and that these equilibria correspond to the steady states of system (HU). It is of crucial 
importance for our analysis that £ is a submanifold of A4 of dimension N + 2, cf. Lemma 2.1 
in [9], Proposition 6.4 in [8]. Now we are prepared to state the main theorem of this paper: 

Theorem 4.2 There exist 7 , S > 0 such that, given po G ®b( 0 , 7 ) D M, problem (14.91) admits 
a unique global strong solution p = £ + e, where (£, e) € E(<S) x £. Moreover, po >—t (£, e) £ 
C 1 (B J5 (0,7)n.M,E( ( 5) xf). 


5 Spectral analysis and proof of the main result 

In this section we study properties of the operator L + K : D(L ) C E 0 > E 0 , where Lp := Lp , 
Kp := Kp and 


D(L) := { (p±,p)eE 1 ; 

a+d n p+ + lpj + Ap = 0 on S , 

[a<9 n /x]] =0 on S, 

d n p _ = 0 on dC. } 

We will identify operators and vector spaces with their complexifications without further men¬ 
tioning. 

Lemma 5.1 (i) The spectrum ofL+K consists purely of isolated eigenvalues having eigenspaces 

of finite dimension. 

(ii) The value A = 0 is an eigenvalue of L + K with dim7V"(L + K) = N + 2. 

(iii) All other eigenvalues of L + K are real and negative. 

PROOF: (i) The operator L generates a strongly continuous analytic semigroup on Eq (The¬ 

orem 2.2 in HD- As the operator p ^ s'(0)p|s • n is of order 1 , K is a relatively compact 
perturbation which implies that L + K also generates a strongly continuous analytic semigroup 
on Eq and in particular has a nonempty resolvent set. Since D{L ) is compactly embedded in 
Eo, the statement follows by [7]’ Theorem III. 6.29. 

(ii) We introduce the bilinear form (•, •) by 


((w+, W-, a), (m_|_, U-, d)) := c_ / w+u+ dx + c+ / W-U-dx^c+C- / aAddS. 

Jd + Jd _ J s 

Letting u(a) := s'( 0)er we observe that 

-({L + K)(w+, w~ ,a), (w+, W-,d)) 

= — c_k+ / A w+w+dx — c+K- / Aw-W- dx + c+c-a+ / d n w+(a+d n w+ + [to]) dS 

Jn + dn_ Js 

+ / u(cr) ■ nAd dS 

Js 

= c_k + ||Vu> + ||! 2 (q + ) + c + ac_||Vm;_ ||| 2 (q_) + n + a + c-\\d n w + \\ 2 L 2^ 

+ i( / \e(u(a))\ 2 dx + [ \e(u(a))\ 2 dx) > 0 (5.1) 

Jd + Jd _ 

for all (w+,w-,a) G D{L). Suppose (L + K)(w+,w-,a) = 0 for (w+,w-,cr) G D(L). Then 
(ED implies that w± are constant on D± (thus [w] is constant on S) and that u{a) = 0 in C. 
From elementary properties of A s and the results from Sections 17.1.11 17.1.31 we hence get that 

A f{L + K) = span{(m, 0, -1), (0, m, 1), (0,0, xi),..., (0,0, xat)} =: {ei, ..., ejv+ 2 }- (5.2) 
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(iii) The computation (EB shows that —((L + K)(w + ,W-,a),(w + ,W-,a)) > 0 for all 
(w+,w-,a) £ D(L). Using the fact that f s u(a ) ■ n = 0, we can show in completely the 
same fashion as in the proof of Lemma 3.1 iii) in [9] that ((w+, w _, a), (w+, w~, cr)) > 0 for all 
eigenvectors (w+,W-,a) of L + K with equality only if (to+ , W- , a) £ A f(L + I\). Hence, the 
assertion follows from Lemma o and (i). ■ 

Corollary 5.2 We have 

sup Re(cr(L + K ) \ {0}) < 0. 

Lemma 5.3 The zero eigenvalue of L + K is semisimple, i.e. X = Af(L + K) © 1Z(L + K). 

Proof: As L + K has nonempty resolvent set and D(L) is compactly embedded in Eq, L + K 

(considered as a bounded operator from D(L) to Eq) is Fredholm and has index zero. Hence, 
by Corollarv l7.6l it suffices to show that J\f(L + K) n 1Z(L + K) = {0}. We introduce the linear 
mapping 


$(/+,/-, 0) : = (c+ Js 0ds + f D+ f+ dx ’ 5 - Is 9dS - Id. f~ dx ) 

and verify assumptions i), ii) and iii) of Lemma 17.31 with 

((w+,w-,8 )\(v+,v-,t)) = {(w+,w-,e), (v + ,v-,f)) 

and V = C 2 . Observe that by (15.11) (and the considerations below (15.11) 1 we have that 

((L + K)(w+, W-, a), (w + , w-, a)) = 0 iff (u)+, u)_, a) £ J\f(L + K) 

The divergence theorem and the fact that j s u(9)-ndS = 0 imply that $ vanishes on 1Z(L + K). 
Assume that 2 := (■ w + , W-,a ) £ Af(L + K) (~l A/"($), (w+,w-,a) = (cf. 15.21) . Then, 

as f s Xj dS = 0 (using a± = k±/c±), $>z = 0 means that 


[ (a 2 - ai) dS = - [ 
Js Jl 


a\m 


d+ c + 


and 


/ («2 — ai) dS = + 

Js Jl 

This is equivalent to A(ai,a 2 ) = 0, where 


a 2 m 
d _ c~ 


dS 


dS. 


( -m\D + \/c+ + \S\ -|5| 

V +\S\ m\D_\/~c_-\S\ 

We calculate, using |5| = iV|D+|, [c] = m = N — 1 

det(j4) = -m\S\\D.\ + m|g||£>_l + m\S\\D + \ + m|£> + ||L>_, > ^ 


Nc- 

i.e. a\ = a 2 = 0. Hence, z £ O. 


C- 


c+ 


Cl¬ 


io 








We turn to the proof of Theorem 14.21 Lemmas 15.11 15.31 allow to follow the same strategy as 
in [9], Section 4. In fact, the arguments given there can literally be repeated here if one only 
replaces the operators L , L by the operators L + K, L + K, respectively. However, the proof of 
[9J, Lemma 4.3 has to be modified because Theorem 2.2 in [4] does not apply to the nonlocal 
operator L + K. Nevertheless, the analogous result holds true: 

Lemma 5.4 Let u> > 0. Then 

(w - (L + K),B) G C is (E,W~ 2 /P(D ± ) x lW'- 3 'P(S)} 3 )nC is (Ei,L p (D ± ) x [W^iS)} 3 ). 

In accordance with our general notation, we denote by Wp 2 ^ p (D±) := Bpp^ p (D±) a Besov 
space of negative differentiability order, see [13- 

Proof: From Lemma 4.3 in J5] we know that 

(w -L,B)e C is (E,W~ 2 /P(D ± ) x [W p 1_3 / p (S')] 3 ) n L P (D±) x [^"^(S)] 3 ). (5.3) 

In particular, (u> — L, B) is Fredholm and has index 0. Since K is a compact perturbation, the 
same is true for the operator (u> — (L+K),B). Therefore it suffices to show that ( ui — (L+K),B ) 
is injective. 

We first consider this operator as an element of C(Ei,L p (D±) x [VFp” 1 ^ p (5)] 3 ). Then, 
injectivity is a direct consequence of Lemma 15. II 
To prove the remaining part, assume 

{u-{L + K),B)n = 0 , p = (p±,p) € E, 


or equivalently 

(w - L, B)p = (0,0, s'(0)p| s • n, 0,0). (5.4) 

Recall that s'(0)p is defined by the BVP (14.81) . Applying Theorem 17.11 to this problem and 
using p G Wp~ 3 P (S), p > 2, we find that the right side of (15.41) is in L p (D ±) x [IFp~ 1 ^ p (S')] 3 , 
so m yields p G E\. The result follows again by Lemma T5. II ■ 

6 Conclusion 

Our analysis crucially relies on the fact that the problem under consideration belongs to the 
class of parabolic evolutions, in the general sense that the semigroup of operators (on ap¬ 
propriate function spaces) arising as solution of the linearized evolution problem is analytic. 
Corresponding maximal regularity results allow the treatment of the nonlinearities introduced 
by the transformation to a fixed domain. As typical for the techniques used here, they provide 
smooth solutions but are (in absence of further structural information) restricted to “pertur¬ 
bative” results, producing either short-time solutions (as in [10J for the present problem) or 
long-time solutions near equilibria or periodic solutions. 

The present paper shows that these techniques are strong and versatile enough to treat 
relatively complex models in which coupled evolutions in two phases and on their interface as 
well as additional elliptic systems occur. On a technical level, this is reflected in the fact that we 
use products of spaces of functions with different domains of definitions and a solution operator 
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for the Stokes equations. In a sense, using this solution operator allows to treat the present 
problem as a perturbed version of the problem without flow, with the perturbation being “of 
lower order.” 

The convergence result may be viewed as an application of a suitably generalized principle 
of linearized stability to a nonlinear parabolic problem, which is also well established by now. 
Discussing the spectrum of the linearization at an equilibrium provides additional structural 
information to conclude that a solution starting close to the manifold of equilibria is actually 
global and converges to this manifold at an exponential rate. 

In this respect, it remains an open and interesting question whether, and how, structural 
properties like parabolicity and stability of equilibria can be concluded already from properties 
of the initial ingredients of the variational model, and not only from the resulting moving 
boundary problem. 


7 Appendix 

7.1 Two-phase Stokes equations 

Let C C R jV be the set defined in the introduction. In this section we denote by 17+ a bounded 
simply connected open set with boundary <917+ of class C°° such that 17+ C C and define 
17- := C \ 17+. Moreover, n denotes the outward unit normal field of 917+. If no confusion 
seems likely, the symbol d n stands for both the directional derivative w.r.t. n and w.r.t. the 
outer unit normal field of dC. We are interested in the two-phase Stokes system 


v±Au± + Vp± 

= f± 

in 17+, 

—div u± 

= 9± 

in 17+, 

[r(u,p)]n 

= h 

on <917+ 

M 

= l 

on <917+ 

U- 

= 0 

on dC 


and consider first the question of unique solvability of the simplified problem 


u±Au± + Vp± 

= f± 

in 17+, 

—div u± 

= 9± 

in 17+, 

[r(w,p)]n 

= h 

on <917+, 

M 

= 0 

on <917+, 

U- 

= 0 

on dC. 


(7.1) 


(7.2) 


7.1.1 Weak solutions 

Let H := Hq(C, R n ) = {w G W 2 (C, 1^); w=0on dC }, Q := L 2 {C), 


a:= H xH 


(u,p)>-+ / v(e{u) : e((p)), 
Jc 


v = v± in 17+, and 


b := Q x H —5- R, (q, ip) i —> — / gdiv<p. 

Jc 
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A weak solution of the system \T?2\ ((/, g) : C —>• R^ -1-1 , h : <912+ — ► R N ) is a pair (u, [g]) G 
H xQ/ ~ c (/ ~ c 5 :<=>/ = g + const) that satisfies 

a(n, <p) + b(q, tp) = / hip + ftp + Vgp for all tp G H (g G [g]) 

9c 9c 

as well as 

bftp, «) = g for all if) G Q. 

Due to Korn’s inequality, the bilinear form a is coercive on H . Moreover, the bilinear form b 
induces a linear operator 

B : H —>• Q', Bu(i/j) := b(xp,u). 

Identifying Q with its dual by means of the Riesz isomorphism the range of B is the set 
{r G L 2 {C ); f c r = 0}. Since this is a closed subset of L 2 (C), classical results (cf. [T], Section 
II. 1) imply that there is a unique weak solution of (17.21) for every (f,g,h) G L 2 (C, R N ) x 
77 1 (C, R) x L 2 (dfl + ,R N ) provided f c g = 0. 


7.1.2 The Lopatinskii-Shapiro condition 


We want to show that the two-phase Stokes system OD satisfies the Lopatinskii-Shapiro con¬ 
dition. W.l.o.g. we restrict ourselves to the halfspace situation, i.e. the case 12+ := R ^ -1 x RA 
Reflecting u- in system EH) to the upper half space, the operator on the left hand side of 
system EH can be expressed by the (2N + 2) x (2N + 2) and the (2 N) x (2 N + 2) matrices 
of operators 


A{dd N ) 


' — z; + A V 


v : 

0 

0 

0 0 | 
—v- A V 1 

, W,- 

••^v) : =( B i + -o N ~ B r e o N ) 

V 0 

0 

V T o / 




both acting on vectors (it+,p+, u_,p_) T (u± := (u+,..., it+)). Here, we used the N xiV-matrices 
of operators 


/a 

o\ 


( ±<9/v 

0 

01 ^ 

A:= • 


-H 

a. 

II 

-H 

m 



i 

Vo 

. 

A/ 


±9iv 

dlV-1 / 
zL25jv / 


and the notation V := (<9i,... ,8n- i, — 9at) t , eN := (0,..., 0,1) G R N . The operator A repre¬ 
sents a Douglis-Nirenberg elliptic system (cf. [3]) with DN-numbers 


SI = ... = SN = SjV+2 — = S2AT+1 = tl = ... = tjV = tN +2 = ••• = t 2 N +1 = 1, 


SN+1 = S2AT+2 = tN+1 = t-2N+2 = 0 

and it coincides with its principal part (note that Sj+tj = 4 N = ord(A)). The characteristic 
polynomial is v+ _1 k^ _1 (|£| 2 + \ 2 ) 2N , where £ = (£ 1 , ..., £jv- i) G R^ -1 and A G R. 

We have to determine a 2AC-dimensional space A4° of exponentially decaying solutions to 
the initial value problem A(i£i ,..., *£jv-i, c 2 t)(u+,p+, u_,p_) = 0, i.e. (letting v± := u+) 

( v ± {\f\ 2 - d^)u° ± , = —ifjP± j = 1, 

{ v±m 2 -d 2 )v ± , = T d P± (7.3) 

l = ±idv±. 
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This system can be solved to the result 


= a£e 1 ^ 1 * 

v± = ]fy(a±|£) 

P± = 0 


j = 1 ,N — 1, o4 G R, 

e _ls|t , a± := (a±, a^' 1 ) 


and 


w± = /M(]fp - 7§f-)e l4|t , j = l,...,iV-1, 

P± = 2 /J ± iy ± e _l4|t , 

/3± G M. Next we show that the problem 

B(i£ i, ...,zYv-i,<9 t )(u + ,p + ,zz_,p_)| t=0 = 0 


(7.4) 


(7.5) 


possesses in M° only the trivial solution. Writing simply (zi+,..., u+ 1 , z> + ) instead of (zzj_,..., u± 


N- 


(u)j_, ...,u± 1 ,u±), we assume that 


«>•■■>< *,«+) - (zA,...,uA 1 ,v_) (0) = 0. 


(7.6) 


This implies that 


J rJ i 


lei 

Multiplication with £ yields 

If, additionally, 


+ 3+-/3_) = 0, (j = 1,...,JV — 1), (a+|0 =(7-7) 


v+ 


/ dtu\ 


9 ‘“+ 
\ d t t; + 




*£iv-iu+ 

dtv+ 


{*± |£) = ^(£+-/?-)• 


/ -Stub 


*£iu- 


LV 


-9t«- 


we find that 




Y -“+l€l > 

+ ( 

A -®(“+k)> 

+ { 

17 «_|C| > 

+ ( 

■U7’ 

1 

7 


( 


-d t v- 


P+ J 


0 

, P+-P- 


(7.8) 

( 0 ) = 0 , 

(7.9) 


+ ( / 3 + ) 


-*(a+IC) J 

t( a -\£)/\t\ \ _ 

-*(«-!£) J 

- 

= 0 . 

By multiplying the first line in (17.101) with £ and by using (17.81) this gives the linear system 


(7.10) 


w -=(::i:z -K7+i-> 


)( 


/?+ 


= o, 


(7.11) 


1 w±)+ 
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which possesses only the trivial solution /3+ = /3_ = 0 since det(M) = — 2(is + + is_) 2 < 0 . 
Hence, by (17.71) . (17.81) a + = a_ and (a + |£) = (a-|£) = 0. Using this, the first line in (17.101) 
reduces to 

(z/+ +i/_)|£|a + = 0 

and hence also a+ = a_ = 0. Therefore, the Lopatinskii-Shapiro condition is satisfied. 

7.1.3 Regularity 

We are now interested in strong/classical solutions of the system (17.11) under the necessary 
solvability demand 



A v+ ' V _(u+,p+,u-,p-) := 13 - 

(y + Au + — Vp + , V- Au- — Vp_, — divu+, — divu_, [r(u,p)]n, |u]|). ^ ' 

From Theorem 9.32 in [3] and Section 17. 1.21 we know that the operators A„ +! „_, considered as 
bounded operator between appropriate function spaces (see Theorem 17. II below h are Fredholm 
for all positive v +, V -. In order to calculate their index, we first consider the case that = 
z/_ =: v > 0 and determine the range of A^,, (for the sake of brevity we refrain from stating 
regularities as they can be easily added by means of classical elliptic theory and the results 
from [2], Section 3.3, 3.5): let l+,l~ satisfy 

• l + — l_ = l on (9H + ; 


• fn± 9± — T fg n+ 4 • n. 


Takine into account (17.121). one nossible choice is L i_ 

:= l + l- |an + , where l- := 

VT and L solves 

-A L = 

g- 

in fl_, 


d n L = 

fn 9- 

on 

(7.14) 

|3fi+| 

d n L = 

0 

on dC. 


Further, let w± := VW±, where W± solve 




-A W+ = 
d n W+ = 

g+ 

1+ ■ n 

in H_|_, 
on dFt+, 

(7.15) 

-AW- = 

g- 

in fl_, 


d n W- = 

l- ■ n 

on <9f2+, 

(7.16) 

d n W- = 

0 

on dC. 


Note that —divu>± = g± in fl ±. We extend /_ and A W- to \ in such 

vanish outside some open ball containing C and consider the problems 

a way that they 

—vAv+ + V( 7 + = 

/+ + vA w+ in f2 + , 


div v + = 

0 

in 

(7.17) 

v+ = 

: 4- 

w + on 



15 









(7.18) 


—vAv- + Vq- = f~ + vAw- 
div v- = 0 

V- = l- — W- 


vN 


m 
in 
on dfl 


\ fi+i 

N \n 


+i 


It follows from Section 3.5 in [2] that the problems (17.171) and (17.181) possess classical solutions 
(since w± ■ n = l± ■ n on <9S1±). Moreover, f gc v~ ■ ngc = 0. Next we are interested in the 
system 


—vAu- 


+ Vp+ 

= 0 

in fi + , 

+ Vp_ 

= 0 

in R N \n+, 

div u -f 

= 0 

in fi + , 

divu_ 

= 0 

in r n \ n + , 

[u,p)Jn 

= h - {t(w + v,q)jn 

on 

M 

= 0 

on (9H + , 


(7.19) 


where v, w : U —> 1^, q : fi + U —> R. are defined in the obvious way. The single layer 

potential with density ^ and w.r.t. the constant viscosity v > 0 is given by 


V i X >'l>) ■= 2^ SA ( n -2)\x-y\»~* + ^ \x-v^ )i>{y) My )i 

1 


\x-y\ r 


Q( x ,fp) ■= zb d<j{ y ). 

As it can be seen from the results in [2], Chapter 3, the restrictions (u±,p±) of 


(7.20) 


(V{-,h- lT(w + v,q)}n),Q(-,h- fr(w + v, g)]n)) 

to 0_)_ and R N \ f2 + , respectively, solve (17.191) in a classical sense, provided h — [ t(w + v, g)]n 
is continuous (observe that precise regularity properties of (u±,p±) can be obtained from the 
fact that the mapping if) k* V(-,i/j)\dn + is a pseudodifferential operator of order —1 as well as 
regularity theory for the Stokes-Dirichlet problem, cf. Section 3.3, 3.5 in i). 

Since u+ is divergence free, it follows that f dn+ u+ ■ n = 0, and [u]| = 0 on implies 
f dn+ U—-n = 0. Hence, since also u_ is divergence free, it follows that J gc U- ■ ngc = 0. Thus, 

fac( u - +l ’A ' n 9C = 0 . 

Let ($, P) a (smooth across <9fi+) solution of the following Dirichlet problem for the Stokes 
equations (which exists since f gc (u- + u_) • ngc = 0, cf. [2], Chapter 3): 


-i/A<L + VP = 0 in C, 

div <I> = 0 in C, 

d> = — (u- + V- + W-) on dC. 


(7.21) 


Summarizing, the pair 

(w± + v± + u± + $, q± + p± + P) 

is easily seen to solve (17.11) (with v + = = v) in a classical sense. Therefore, the necessary 

solvability demand (17.121) is also sufficient. Hence, the range of „ is of codimension 1 . Since 
we know from Theorem 1 7.1.1 1 that the kernel of is one dimensional, this operator has index 
0. Consequently (y > 0 was arbitrary), by homotopic stability of the index, all members of the 
family {A v+ Ji-t)v + +tv_> t € [0,1]} have index 0, in particular A„ +j „_. Since also this operator 
has a one dimensional kernel, the following theorem 1 7. II can be deduced from the general theory 
of elliptic boundary value problems (cf. Section 4 in [17] . Theorem 9.32 in 0). In order to 
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economize notation we introduce the quotient spaces T := T / ~ c , where T G {W“, G a , c a } 
and ~ c is the equivalence relation introduced in Section 17.1.11 Here, C a stands for the usual 
Holder space and c“ denotes the little Holder space, that is the closure of the smooth functions 
in C a . 

Theorem 7.1 Let (r, /3,p, k) G (0,1] x (0,1) x [1,oo) x(MU {0}) satisfy k + r > 1/p, r/ 1/p 
and let C G {c,C}. Suppose that ( f±,g±,h,l) of class 

i) W := [TTp +r_ 1 (f2±)] Ar x W k+r (fl ± ) x [Wp +r ^ 1/p {dfl+)} N x [Wp +r+1 ~ 1/p (dfl + )} N 

or 

a) c ■- x c k+1+ p(£i±) x [c k+l +v(dn+)] N x [c k+2+f> (dn+)] N 

satisfies f n g+ + g~ = — l-n. In both cases problem (17.11) possesses a solution ( u±,p ±) 
which is unique up to an additive constant for p±. This constant is the same in both phases. 
In case i) ( u±,p ±) belongs to the class [H / p +:r+ 1 (H±)] Ar x W k+r (12±) and satisfies the a priori 
estimate 

ll( u ±> [l'±])ll[Wp fc+, ' + 1 (n ± )] Jv xW' p fe + r (n ± ) — 7ll(/±)5±>^) Ollw (7.22) 

with a positive constant 7 independent of (f±,g±,h 1 l). In case ii) (u±,p±) belongs to the class 
[C k+2+ P(Ll±)] N x C fc+1 +^(H±) and satisfies the estimate 

ll(«±,b±])ll [c fc + 2 +^(n ± )] N x c fc +i+/ 3 (n±) <7ll(/±,5±,^0llc (7.23) 

with a positive constant 7 independent of (f±,g±,h,l). The symbol [p-j_] stands for the class 

{p± + Ci C G R}. 

For a given domain C C with smooth boundary we denote by Aq ± the operator defined in 
(17.131) . Using this notation (and those defined in Section [4j, we are able to work out that the 
solution operator s(p) of the problem (14.31) depends smoothly on p. 

Corollary 7.2 There is a neighbourhood U of 0 in Wp 3 ^ P (S) such that s(p) is well defined 
for all p G U D Ad and s G C°°(U D Ad, {Wp~ 2 ^ p (D±)] N ). Moreover, if 5 > 0 is given, then 
there is a neighbourhood V of 0 in E(<5) such that s given by 

s((g±,cr))(t) := (0, 0, s(cr(t))|s • n s ) 

is well defined for all ( p±,cr ) G V and s G C°°{y, L p (K + , Eq)). 

Proof: Let A (p) := 0*A(Q p ) ± 0*. Carrying out the transformation of the differential opera¬ 

tors involved in (14.31) (cf. the proofs of Theorem 3.1 in [TU] . Lemma 4.2 in 0 ) we obtain 

s(p)=PA(p)- 1 (0,0,i7(p)n(p),0), 

where P{u+,p+, U-,p-) := u and A, H and n depend all smoothly on p. Since A(0) is precisely 
the operator considered in Theorem 17.11 (with II± := D±), the first assertion is obtained by 
standard perturbation arguments for isomorphisms. 

Observe that due to the embedding (14.71) we may assume that er[R + ] C Ad D Wp~ 3 ^ p (S) 
provided |[(/i±, cr) |]e( 5 ) is small enough. Hence the second assertion is a consequence of mapping 
and smoothness properties of concerning Nemytskij operators, cf. m- ■ 
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7.2 A few functional analytic tools 

Lemma 7.3 Let H , V be complex vector spaces and A : D{A ) C 77 —> 77 be a linear operator. 
Let (-I-) be a sesquilinear form on 77(A) x 77(A), $ e £(77, V) and 

O := {z £ 77(A); ( z\h ) = 0 for all h £ 77(A)}. 

Assume that 


i) (Au\u) = 0 iff u £ Af(A); 

ii) 77 (A) c 7V($); 

in) M($)nAT(A) C O; 

Then 77(A) nTV(A) = {0}. 

Proof: Let 2 : € 77(A) nA/"(A), z = Ay for some y £ D(A). Then ( because of ii), iii) ) z £ O 

and 

0 = (z\y) = (Ay\y), 

and hence, because of i), y £ Af(A) i.e. z — Ay = 0. ■ 

Lemma 7.4 Let H be a complex vector space, A : D(A) C 77 —>• H and let (-|-) be a sesquilinear 
form on D(A) x 77(A). Assume additionally that 

i) {Au\u) £ [0, 00 ) for all u £ 77(A); 

ii) («|u) £ [0, 00 ) for all eigenvectors u of A; 

iii) if (u\u) = 0 for some eigenvector u of A, then u £ AT (A). 

Then all eigenvalues of A are nonnegative. 

Proof: Let A ^ 0 be a (possibly complex) eigenvalue of A. Then, if u £ 77(A) is a cor¬ 

responding eigenvector, we have by ii) and iii) that (u\u) > 0. The assertion follows from 
A(it|it) = (Aw|it) = (Au\u) >0. ■ 


Lemma 7.5 Let X be a Banach space, N,R linear subspaces with the properties that TV D 77 = 
{0}, dim(N) = codim(R) = M < 00 and that R is closed. Then the quotient map Q : X —> X/R, 
y n- y+R induces a topological isomorphism from N onto X/R and it vanishes on R. Moreover, 
X = N © R algebraically and topologically. 

Proof: It straightforward to check that Q|jv is a topological isomorphism onto the finite 

dimensional Banach space X/R and that Q vanishes on R. Moreover, P := (Q|at) _1 o Q is a 
continuous projection of X onto TV, hence 

X = 1 V © Af(P), dim(J K/Af{P)) = M. 

Since R C Af(P) and dim(X/7?) = M, we find dim(Af(P)/R) = dim(X/R) — dim(X/Af (P)) = 0 
and thus Af(P) = R. ■ 

Since the range of a Fredholm operator is always closed, we get the following 

Corollary 7.6 If X,Y are Banach spaces such that Y C X, and if F £ C{Y. X) is a Fredholm 
operator of index 0 that satisfies 77(F) C\Af(F) = {0}, then X = Af(F) ©77(F) algebraically 
and topologically. Moreover, P := {Q\j£(f))~ 1 o Q is a continuous projection of X onto Af(F). 
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